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[i] The high time resolution observations obtained by the WAVES experiment of the 
STEREO spacecraft in solar type III radio bursts show that Langmuir waves often occur as 
intense localized wave packets. These wave packets are characterized by short durations of 
only a few ms and peak intensities, which well exceed the supersonic modulational 
instability (MI) thresholds. These timescales and peak intensities satisfy the criterion of the 
solitons collapsed to spatial scales of a few hundred Debye lengths. The spectra of these 
wave packets consist of primary spectral peaks corresponding to beam-resonant Langmuir 
waves, two or more sidebands corresponding to down-shifted and up-shifted daughter 
Langmuir waves, and low frequency enhancements below a few hundred Hz corresponding 
to daughter ion sound waves. The frequencies and wave numbers of these spectral 
components satisfy the resonance conditions of the modulational instability (MI). 

Moreover, the tricoherences, computed using trispectral analysis techniques show that 
these spectral components are coupled to each other with a high degree of coherency as 
expected of the MI type of four wave interactions. The high intensities, short scale lengths, 
sideband spectral structures and low frequency spectral enhancements and, high levels of 
tricoherences amongst the spectral components of these wave packets provide 
unambiguous evidence for the supersonic MI and related strong turbulence processes in 
type III radio bursts. The implication of these observations include: (1) the MI and 
related strong turbulence processes often occur in type III source regions, (2) the strong 
turbulence processes probably play very important roles in beam stabilization as well as 
conversion of Langmuir waves into escaping radiation at the fundamental and second 
harmonic of the electron plasma frequency, f pe , and (3) the Langmuir collapse probably 
follows the route of MI in type III radio bursts. 

Citation: Thejappa, G., R. J. MacDowall, and M. Bergamo (2012), In situ detection of strong Langmuir turbulence processes in 
solar type III radio bursts, J. Geophys. Res., 117 , A081 1 1, doi:10.1029/2012JA017695. 


1. Introduction 

[ 2 ] This paper describes the strong Langmuir turbulence 
processes, namely, the supersonic modulational instability 
(MI) and collapsing Langmuir solitons observed in the 
source regions of solar type III radio bursts. These obser- 
vations are from the WAVES experiment of the STEREO 
spacecraft [Bougeret et al., 2008], which is designed to 
detect in situ waves as well as remote radio emissions over a 
wide frequency range. Ginzburg and Zheleznyakov [1958] 
were the first to suggest that (1) the electron beams accel- 
erated during solar flares generate Langmuir waves at the 
local electron plasma frequency, f pe [Hz\ = 9 n l J 2 by the 
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bump-on-tail instability [Bohm and Gross, 1949] ( n e is the 
electron density in m~ 3 ), and (2) these Langmuir waves are 
subsequently converted by some nonlinear processes into 
the observed type III radio emissions at the fundamental and 
second harmonic of the local electron plasma frequency, f pe . 
Although this hypothesis has been confirmed by the in situ 
detection of type III burst associated electron beams [Lin, 
1970; Lin et al., 1973, 1981; Ergun et al., 1998], and 
Langmuir waves [Gurnett and Anderson, 1976, 1977; Lin 
et al., 1986; Gurnett et al., 1993; Thejappa et al., 1993a, 
1993b, 1999; Hospodarskv and Gurnett, 1995; Thejappa 
and MacDowall, 1998; Henri et al., 2009], several issues 
remain unresolved. 

[ 3 ] Sturrock [1964] was the first to pose the following 
question: How does the electron beam preserve the bump- 
on-tail distribution over 1 AU and beyond against the 
“quasi-linear relaxation” or “plateauing”, that is known to 
deplete all the streaming energy within 1 00 km or less? Here 
we should note that Sturrock considered a special case of an 
infinitely extended homogeneous plasma, and neglected the 
effects of re-acceleration of beam particles by Langmuir 
waves. Several authors [Zaitsev et al., 1972, 1974; Takakura 
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Figure 1 . Schematic picture of the modulational instability 
(Ml). Two beam excited Langmuir waves L are scattered 
into down-shifted (Stokes) L D and up-shifted (anti-Stokes) 
L v sidebands through the exchange of ion sound waves 
(IS). The frequencies and wave numbers of the beam excited 
Langmuir, and daughter ion sound, Stokes, and anti-Stokes 
modes are (f L , k L ), (O, q), (f L - O, k L ~ q), and (f L + Q, 
k L + q), respectively. 

and Shibahashi, 1976; Magelssen and Smith, 1977; Escande 
and de Genouillac, 1978] have treated the wave particle 
interactions self-consistently all along the path of the inho- 
mogeneous beam, and claimed that the continuous emission 
and reabsorption of waves enable the beam to survive over 
the distances of order of ~1 AU and more. However, the 
Langmuir wave intensities predicted by the inhomogeneous 
quasi-linear theories are also subject to wave-wave interactions. 

[ 4 ] In the context of homogeneous beam-plasma interac- 
tion, Kaplan and Tsytovich [1968] suggested that the 
induced scattering of Langmuir waves off ion clouds, which 
acts as the electrostatic decay of Langmuir waves into 
backscattered Langmuir waves and forward-going ion sound 
waves when T e > T t (T e and 7] are the electron and ion 
temperatures, respectively) can stabilize the beam; these 
authors considered only the mono energetic beam and iso- 
tropic distribution of waves. Although some signatures of 
the electrostatic decay are detected in type III source regions 
[Lin et al., 1986; Gurnett et al., 1993; Thejappa et al, 1995, 
2003; Thejappa and MacDowall, 1998; Hospodarsky and 
Gurnett, 1995; Henri et al., 2009], if the assumption of 


mono energetic beam and isotropic distribution of waves is 
relaxed, the net scattering of waves out of resonance with the 
beam by the induced scattering becomes too slow to prevent 
quasi-linear relaxation [Zheleznyakov and Zaitsev, 1970; 
Hevvaerts and de Genouillac, 1974]. 

[5] The supersonic MI, which is popularly known as the 
oscillating two stream instability (OTSI) [Papadopoulos 
et al, 1974; Smith et al, 1979; Goldstein et al., 1979] and, 
soliton formation and Langmuir collapse [Zakharov, 1972; 
Nicholson et al., 1978; Goldman, 1983] are suggested as the 
most effective beam stabilization mechanisms. These strong 
Langmuir turbulence processes are the focus of several the- 
oretical studies [Zakharov, 1972; Goldman, 1 984; Robinson, 
1997], computer simulations [Bardwell and Goldman, 
1976; Nicholson et al., 1978; Degtyarev et al., 1980; Doolen 
et al., 1985; Russell et al., 1988] and experimental investi- 
gations [Whelan and Stenzel, 1981; Michel et al., 1982; 
Wong and Cheung, 1984; Cheung and Wong, 1985a, 1985b; 
McFarland and Wong, 1997, 2001; Vyacheslavov et al., 
2002] because of the fundamental nature of the problem 
[Zakharov, 1972; Krasnoselskikh and Sotnikov, 1977], and 
practical importance in fields, such as the laser heating of 
pellets [Montgomery et al., 1999], ionospheric modification 
experiments [DuBois et al., 1993], and several astrophysical 
applications [Galeev and Krasnoselskikh, 1976, 1978; 
Gurnett et al., 1981; Pelletier et al., 1988; Balikhin et al., 
1989; Kellogg et al., 1999], 

[6] As shown in Figure 1, MI excites a low frequency ion 
density perturbation of frequency and wave number (Q, q), 
which can be either a quasi-mode or a freely propagating ion 
sound wave. This ion sound wave can beat with two of the 
beam-resonant Langmuir waves with frequencies and wave 
numbers {f L , k L ) and produce high frequency sidebands in 
the non resonant region with frequencies and wave numbers 
(f L — Q, k L — q) (Stokes mode) and (f L + O, k L + q) (anti- 
Stokes mode). The nonlinear effect responsible for MI is the 
density cavity dug by the action of the ponder motive force 
of the Langmuir wave electric field, which in turn traps the 
waves thus deepening it still further. Usually, the Langmuir 
wave bunch localized in the corresponding density cavity is 
referred to as the Langmuir soliton. The spatial collapse, on 
the other hand, occurs due to intensification of the Langmuir 
wave packet in the self generated density cavity [Zakharov, 
1972], 

[7] Gurnett et al. [1981] were the first to report the Voy- 
ager observations of the possible evidence for MI and spatial 
collapse of Langmuir waves in the Jovian foreshock region. 
However, in this case, the daughter ion sound waves were not 
detected, and the intensities of Langmuir waves were well 
below the strong turbulence thresholds. The Ulysses obser- 
vations of Langmuir waves in the form of intense millisecond 
spikes were interpreted in terms of Langmuir collapse 
[Kellogg et al., 1992; Thejappa et al., 1993b]. However, in 
this case, there were some instrumental uncertainties 
regarding millisecond spikes. We have reported the obser- 
vations of the coexistence of strong and weak Langmuir 
turbulence processes [Thejappa and MacDowall, 1998], 
envelope solitons formed as a result of subsonic MI 
[Thejappa et al., 1999] and ion-sound waves radiated by the 
bumt-out cavitons left behind by the collapse events 
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[Thejappa and MacDowall, 2004] in type III sources. How- 
ever, the envelope soliton does not correspond to the concept 
of physical collapse of Langmuir waves because of its rela- 
tively large spatial scale S » A and low intensity (A is the 
wavelength of the Langmuir wave). Of much greater interest 
from the point of collapse is the large-amplitude Langmuir 
soliton formed in the regime J'j > ( k L \ De ) 2 and ^ > ™‘;, with 
S < A = which corresponds to the supersonic MI [Zakharov 
et al., 1985]. Here j-jr is the normalized energy density of 
Langmuir waves, m e and m l are the electron and ion masses, 
X/),, is the electron Debye length and k/ is the wave number of 
the Langmuir waves. Even after years of comprehensive in 
situ high time resolution measurements of the wave and beam 
parameters [Lin et al, 1986; Kellogg et al., 1992; Gurnett 
et al., 1993; Thejappa et al., 1993a, 1993b, 1996, 1999, 
2003; Hospodarsky and Gurnett, 1995; Thejappa and 
MacDowall, 1998, 2004; Ergun et al., 2008; Henri et al., 
2009], the evidence for the Langmuir collapse has proved 
to be extremely elusive. 

[s] Recently, Thejappa et al. [2012a] reported the high 
time resolution observations of a very intense localized 
Langmuir wave packet in the source region of a type III 
radio burst from the improved Time Domain Sampler (TDS) 
of the STEREO/WAVES experiment (improved over that of 
all similar high time resolution receivers flown in earlier 
spacecraft in precision, linearity, sample length and rate 
[Kellogg et al., 2009]). The characteristics of this wave 
packet include: (1) the normalized peak energy density 
which is well above the threshold for the supersonic MI, 
(2) timescale of only ~3 ms, and (3) spectrum containing the 
resonant peak at the local electron plasma frequency, f pe , 
sidebands at frequencies slightly lower than f pe (Stokes 
mode) as well as slightly higher than f pe (anti-Stokes mode), 
and low frequency enhancement below a few hundred Hz, 
which satisfy the frequency and wave number matching 
rules of MI. The spatial scale determined using the observed 
timescale of the wave packet is less than the spatial scale 
derived using the normalized peak intensity, which indicates 
that the wave packet is probably the Langmuir soliton col- 
lapsed to a few hundred Debye lengths. These observations 
were interpreted in terms of the oscillating two stream 
instability (OTS1), which is nothing but the supersonic Ml 
and collapsing Langmuir soliton. In a separate study, 
Thejappa et al. [2012b] applied the trispectral analysis 
technique and showed that the spectral components of this 
wave packet are coupled to each other with a high degree of 
phase coherency (high tricoherence). 

[ 9 ] In this study, we present a detailed analysis of high 
time resolution observations of Langmuir wave packets 
associated with four other type III events. These TDS events 
represent some of the most intense localized Langmuir wave 
packets ever detected in type III bursts. The characteristics of 
these wave packets indicate that Ml and associated strong 
turbulence processes occur in type III radio bursts. The 
implications of this study include (1) the strong turbulence 
processes probably play important roles in beam stabiliza- 
tion as well as conversion of Langmuir waves into escaping 
fundamental and harmonic radio emissions, and (2) the 
Langmuir collapse in type III sources probably occurs via 


MI as suggested by Zakharov [1972]. In section 2, we 
present the observations, in section 3, we present the dis- 
cussion, and in section 4, we present the conclusions. 

2. Observations 

[ 10 ] The STEREO/WAVES experiment [Bougeret et al., 
2008] is uniquely suited to study the nonlinear processes 
associated with solar type III radio bursts. It can detect a 
variety of in situ waves as well as remote radio emissions 
over a wide frequency range. Data from three of the 
SWAVES instruments are pertinent to this study: the high 
frequency receiver (HFR), low frequency receiver (LFR) 
and time domain sampler (TDS). 

[ 11 ] 1 . High frequency receiver (HFR)- The HFR measures 
the remote electric field signals over the frequency range 
125 kHz - 16.025 MHz with 50 kHz increments. The total 
dynamic range of the HFR is 80 dB and its sensitivity at the 
preamplifier input is 6 nVHz -1/2 . 

[ 12 ] 2. Low frequency receiver (LFR)- The LFR is a direct 
conversion receiver for spectral processing from 2.5 kHz up 
to 160 kHz. The frequency range is divided into three 2-octave 
bands A (2.5-10) kHz, B (10-40) kHz and C (40-160) kHz. 
The LFR covers a broad dynamic range with a constant 
amplitude resolution over the whole range, i.e., it adjusts its 
gain according to the input level. Automatic Gain Control 
(AGC) determines the receiver gain on a logarithmic scale as 
a function of the input level. The LFR analog front end nor- 
malizes the signal amplitude with an allowable gain variation 
of 80 dB. Combined with 12-bit analog to digital converters, 
the LFR provides a total dynamic range of 120 dB. The LFR 
sensitivity is 6 nVHz _1/2 at the preamplifier input. 

[ 13 ] 3. Time domain sampler (TDS)- The TDS is intended 
primarily for the study of Langmuir waves at the local 
electron plasma frequency, f pe . It gathers events on three 
channels simultaneously connected to three orthogonal 
antennas and to a fourth pseudo-dipole channel obtained by 
taking the difference of any two monopoles. The TDS has 
programmable sampling rates with events of maximum 
length of 16 kilo samples, with possible rates of 250,000, 
125,000, 31,250 and 7,812.5 samples per second. The 
maximum event duration can be as long as 2 s. The use of a 
16 bit linear A/D converter provides sufficient dynamic 
range with better fidelity. The long durations of the TDS 
events enable us to determine the frequency components of 
the wave packets accurately. The gain is linear in the fre- 
quency range of interest, from typical electron cyclotron 
frequency, f ce ~ 100 Hz to typical electron plasma frequency 
fpe ~ 30 kHz. 

[ 14 ] In the data, we have identified four type III events, for 
which the radio emissions drift from high frequencies to f pe or 
2 f pe , where f pe is the local electron plasma frequency. The 
Langmuir wave packets captured by the TDS during these 
events form the observational basis for the current study. 
These type III events occurred on 201 1 December 19 (event 1), 
2012 February 12 (event II), 2009 July 18 (event III) and 
2006 December 5 (event IV). We assign a typical value of 
~10 5 K for the electron temperature T e for all four events. 
Here we note that owing to an unfortunate malfunction of the 
electron experiment, STEREO does not measure the electron 
temperature, T e . For the solar wind speed (y„ v ), we use the 
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values of 450 km s , 500 km s _1 and 500 km s _1 for events 
I, II and III, respectively. These values are from the STEREO/ 
PLASTIC experiment [Galvin et al., 2008]. Since the data 
of solar wind are not available from the STEREO spacecraft 
for 2006 December 5 event (event IV), we use a value v sw = 
315 km s -1 obtained by the WIND/S WE experiment 
[Ogilvie et al., 1995]; the STEREO A, B and WIND were 
very close to each other during this event. We have 
deduced the electron densities, n e of the solar wind as 
~4.5 x 10 6 m -3 , ~2.4 x 10 6 m“ 3 ,~5x 10 6 m“ 3 and-1.2 x 
10 7 m“ 3 from the frequencies of the dominant spectral peaks 
of the waveforms of ~ 1 9.12 kHz, ~14 kHz, ~20. 1 kHz 
and ^3 1.6 kHz (see, Figures 5, 9, 13, and 17); these fre- 
quencies are assumed to correspond to local electron plasma 
frequencies, f pe [H:\ = 9n\ /2 , ( n e is in m“ 3 ) during events 1 
to IV, respectively. 

[is] The solarsoft solar activity reports indicate that events 
1 to IV are likely associated with active regions of N18W28, 
N29E85, N22W81 and S06E71, respectively. For the elec- 
tron density n e (m~ 3 ), we assume the Radio Astronomy 
Explorer (RAE) density model [Fainberg and Stone, 1971]: 

n e = n 0 r ~ a , (1) 

where « 0 = 5.52 x 10 13 , a = 2.63 and r is the solar altitude 
(in units of Rq). If one assumes that the type 111 bursts are 
excited at 2 f pe by the electron beams traveling along the 
spiral magnetic field lines with velocities (3 (units of velocity 
of light c), the frequency drifts can be expressed as 
[Papagiannis, 1971] 

dJ a P C j-(a+2)/a ( 2 ) 

dt 2(1 — /3cos</>)R 0 (81 x KT 6 ) 17 * 

where c is the velocity of light, <f> is the angle of exciter 
direction to Sun-Spacecraft line, and f is the midpoint of 
frequency interval df. This expression is derived assuming 
that the type III electrons propagate along the Parker’s spiral 
field lines. By fitting the frequency drift curves to the 
dynamic spectra, we can estimate the beam speeds v h as 
~0.3c, ~0.28c, ~0.37c, and ~0.35c for events I to IV, 
respectively. Since, the pitch angle scattering increases the 
path length traveled by electrons by a factor of a = 1 .3 to 1.7 
[Alvarez et al., 1975; Lin et al., 1973; Fokker, 1984], the 
actual beam speeds are probably much greater than these 
estimated speeds. 

[ 16 ] Using these values, we estimate: (1) the Debye 
lengths \ De = 69T’' /2 «J 1/2 as ~ 10.3 m, ^14 m, ~9.8 m and 
~6.2 m, (2) wave numbers of beam-excited Langmuir waves 
k L = ^as x 10“ 3 m _1 , ~1 x 10 _3 m _1 ,~l.l x 10“ 3 
m _1 and ~1.9 x 1(T 3 m _1 , and (3) k L X De as ~1.4 x 10 2 , 
~1.5 x 10“ 2 , ~1.1 x HT 2 , and -1.2 x 10“ 2 , for events 1 
to IV, respectively. 

2.1. December 19, 2011 (Event I) 

[ 17 ] In Figure 2, we present the frequency-time spectro- 
gram illustrating the type III radio burst. The data used in the 
spectrogram are from the HFR as well as LFR of the WAVES 
experiment of the STEREO A spacecraft. The intense emis- 
sion drifting from ~5 MHz to ~28 kHz is the type III burst. 


This is identified as the local event because of its drift from 
very high frequencies all the way to twice the local electron 
plasma frequency f pe . Since type III emission frequencies 
depend on the radially decreasing electron density n e , these 
bursts occur at lower frequencies at later times. The in situ 
wave data acquired during local type III events are exten- 
sively used to study the non-linear plasma processes [Lin 
et al., 1986; Kellogg et al., 1992; Gurnett et al., 1993; 
Thejappa et al., 1993a, 1999, 2003; Hospodarsky and 
Gurnett, 1995; Henri et al., 2009]. 

[is] In association with the type III radio burst, the LFR 
has detected Langmuir waves (non-drifting emissions in the 
19-25 kHz frequency range) generated probably by the solar 
flare electron beam. In Figure 3, we present the frequency- 
time spectrogram of Langmuir waves. As seen from this 
spectrogram, the frequencies of these emissions fall sud- 
denly from ~28 kHz to ~ 1 9 kHz at ~ 13:45 UT, and remain 
at that frequency for almost 35 min, and, then jump suddenly 
to ~25 kHz at ~ 14:20 UT. These abrupt changes in the 
frequencies of Langmuir waves are probably due to some 
density structure in the solar wind. As seen from this spec- 
trogram, at any given instant, the Langmuir waves show 
substantial frequency spreading of ^8 kHz, which is indic- 
ative of nonlinear frequency broadening. The Langmuir 
wave burst corresponding to the TDS event of the present 
study is shown by an arrow in this figure. 

[ 19 ] Time Domain Sampler (TDS) has resolved these 
Langmuir waves as highly localized wave packets. Each of 
these wave packets contains 16384 samples with an acqui- 
sition rate of 250, 000 samples per second (a time step of 4 ps 
for a total duration of 65 ms). In Figure 4, we present one of 
the most intense wave packets detected during this type III 
event. This event is selected because of its sideband spectral 
structure, even though it is slightly weaker than the TDS 
events, which occur later at ^13:20 UT. The peak electric 
field strength E r of this event is 34.6 mVnV 1 and its j - 
power duration r is ~5.7 ms. These values yield (1) the 
normalized peak energy density -pjr = jpp — 8.5 x 10~ 4 for 
E l = 34.6 mVm -1 , n e = 4.5 x 10 % m“ 3 and T e = 10 5 K, and 
(2) the spatial scale of the field structure S ~ 250\ De (using 
the relation S ~ tv sw ) for r = 5.7 ms, v sw = 450 km s -1 and 
\ De ~ 10.3 m. Since the E x , E v and E x _ y signals are weaker 
and show the same general features as the E- signal, we 
analyze only the E : signal. We justify one-dimensional 
treatment by assuming that these Langmuir wavefields are 
probably aligned along the ambient magnetic field. 

[ 20 ] In Figure 5, we show the spectral characteristics of 
this wave packet. Figure 5a shows the total spectrum of the 
waveform from 0 to 65 kHz. The primary spectral peak at 
~19.12 kHz probably corresponds to Langmuir waves 
excited by the beam at the local electron plasma frequency, 
f pe . Figure 5b shows the narrow logarithmic spectrum 
around this primary spectral peak. This spectrum contains an 
intense peak (L) at f pe ~ 19.12 kHz, corresponding to the 
beam excited Langmuir wave, a down-shifted sideband (D) 
at ~19 kHz which is slightly less than f pe and an up-shifted 
sideband (U) at ~19.3 kHz, which is slightly higher than^, e . 
Figure 5c shows linear spectrum below 1.5 kHz. A clear 
enhancement below 200 Hz seen in this low frequency 
spectrum probably corresponds to ion-sound waves. 
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Figure 2. Dynamic spectrum of the wave activity observed by the WAVES experiment of STEREO A 
spacecraft from 1 1:00 to 15:00 UT on 201 1 December 19. The fast drifting emission from ~5 MHz down 
to ~28 kHz is identified as the local type III radio burst, and the non-drifting emissions in the 19-25 kHz 
frequency band are identified as the Langmuir waves. 


2.2. February 11, 2012 (Event II) 

[ 21 ] In Figure 6, we present the dynamic spectrum of the 
second type III event and associated Langmuir waves, 
observed on February 1 1, 2012 by the STEREO B spacecraft. 
The fast drifting type III burst and the non-drifting Langmuir 
wave emissions are clearly seen in this dynamic spectrum. 
In Figure 7, we present the frequency-time spectrogram in a 
narrow frequency range, which shows the clumpy Langmuir 
waves with substantial frequency spreading. In this figure, 
the Langmuir wave burst corresponding to the TDS event of 
the present study is shown by an arrow. In Figure 8, we 
present one of the most intense wave packets detected by the 
TDS during this type III event. This TDS event is charac- 
terized by the peak electric field strength E L of 44.73 mVm“ 1 
and i -power duration r of ^7.9 ms. These values yield: 
(1) nje ~ 2-6 x 10 -3 for E l = 44.7 mVm" 1 , n e = 2.4 x 10 6 
m" 3 and T e = 10 5 K, and (2) S ~ 2%2\ De for the measured 
values of r = 7.9 ms and v sw = 500 km s -1 . In Figure 9, 
we show the spectral characteristics of this wave packet. The 
total spectrum of the waveform from 0 to 65 kHz is shown in 
Figure 9a, where the primary spectral peak at ~ 1 4 kHz cor- 
responds probably to beam excited Langmuir wave at f pe . 
The narrow logarithmic spectrum around this primary peak 
(L) in Figure 9b shows two sidebands, Stokes peak (D) at 
^13.78 kHz ( which is slightly less than f pl 2), and an anti- 
Stokes peak (U) at ~1 4.25 kHz ( which is slightly higher than 
f pe ). In Figure 9c, we present the low frequency spectrum, 


which shows a clear enhancement <250 Hz corresponding to 
ion-sound waves. 

2.3. July 18, 2009 (Event III) 

[ 22 ] In Figure 10, we present the dynamic spectrum of the 
third local type III event and associated Langmuir waves. 
This event was observed by the STEREO A spacecraft on 
July 18, 2009. In Figure 11, we present the frequency-time 
spectrogram of the Langmuir waves, where the arrow shows 
the location of the TDS event of the present study. Here we 
note that even though this TDS event (Figure 12) is slightly 
weaker than some of the later events, it is included in the 
present analysis because of its sideband spectral structure. 
This TDS event is characterized by the peak electric field E L 
of 70.9 mVm 1 and ~ e -power duration r of ~ 1 0 ms. These 
values yield: (1) H ,r '. - 3.2 x 10~ 3 for E L = 70.9 mVm" 1 , 
n e = 5 x 10 6 m -3 and T e = 10 5 K, and (2) 5 - 5\2\ De for the 
measured values of r = 10 ms and v sw = 500 km s^ 1 . 
Figure 13 presents the spectral characteristics of this wave 
packet. Figure 13a shows the total spectrum of the waveform 
from 0 to 65 kHz, where the dominant spectral peak at 
~20. 1 kHz corresponds to Langmuir waves excited by the 
electron beam at f pe . Figure 13b shows the narrow logarith- 
mic spectrum around this intense peak (L), which contains 
two sidebands: a Stokes peak (D) at ^ 19.85 kHz (which is 
slightly less than f pe ), and an anti-Stokes peak (U) at 
~20.32 kHz (which is slightly higher than f pe ). In 
Figure 13c, we present the low frequency spectrum of this 
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Figure 3. Frequency-time spectrogram of the Langmuir wave activity associated with the type III burst 
of December 19, 201 1. These waves are very bursty and show substantial frequency broadening. As seen 
from this spectrogram, the frequencies of these emissions fall suddenly from ~28 kHz to ~ 1 9 kHz at 
~13:45 UT, and remain at that frequency for almost 35 min, and, then jump suddenly to ~25 kHz at 
~ 14:20 UT. These abrupt changes in the frequencies are probably due to some unknown density structure 
in the solar wind. The location of the current TDS event is shown by the arrow. 


STEREO A - 19 December, 201 1 



wave packet, which shows a clear enhancement <250 Hz 
corresponding to ion-sound waves. 

2.4. December 5, 2006 (Event IV) 

[23] In Figure 14, we present the dynamic spectrum of the 
fourth local type III event and associated Langmuir waves, 
observed by the STEREO B spacecraft. Here the fast drifting 
emission from ~ 1 0:00 to ~ 12:30 UT is the type III burst, 
and the non-drifting emissions in the 28-32 kHz frequency 
range are the Langmuir waves. In Figure 15, we present the 
frequency-time spectrogram in a narrow frequency band, 
which clearly shows that the Langmuir waves are very 
bursty and show substantial frequency spreading ~10 kHz. 
The arrow shows the Langmuir wave burst corresponding to 
the TDS event of the present study. 

[24] In Figure 16, we present one of the most intense 

waveforms detected during this type III event, which con- 
tains several modulated structures. The most intense struc- 
ture is characterized by El = 69.7 mVm 1 and r ~ 7 ms. 
These measured parameters yield: (1) — 1.3 x 10“ 3 for 

E l = 69.7 mVm" 1 , n e = 1.23 x 10 7 m J ' j and T e = 10 5 K, and 
(2) S ~ 355X De for the measured values of r = 7 ms and v sw 
= 315 km s^ 1 . In Figure 17, we show the spectral char- 
acteristics of this wave packet, where Figure 17a shows the 


total spectrum from 0 to 65 kHz. The dominant spectral peak 
at ~31.6 kHz in Figure 17a corresponds to Langmuir waves 
excited by the electron beam at the local electron plasma 
frequency,^, which corresponds to n e of 1.23 x 10 7 m -3 . 
The narrow logarithmic spectrum around this dominant 
spectral peak (L) as shown in Figure 17b contains several 
secondary peaks. The peaks at frequencies slightly less than 
fpe (Stokes modes) are denoted as D x and D 2 , and peaks at 
frequencies slightly higher than f pe (anti-Stokes modes) are 
denoted as U\, U 2 and 63 . The low frequency spectrum as 
seen in Figure 17c clearly shows the ion acoustic wave 
associated enhancement <200 Hz. 

3. Discussion 

[25] In Table 1 , we present the summary of observations. 
The normalized peak energy density j- jr = controls the 

non-linear beam-plasma interactions. If 

— — > ( k L \ De ) 2 , (3) 

YleT-e 

then the dominant nonlinear beam-plasma interactions are 
the modulational instability [ Papadopoulos et al., 1974] and 
the fully developed soliton formation and collapse 
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Figure 4. One of the most intense Langmuir wave packets captured by the Time Domain Sampler (TDS) 
during type III event of Figure 2. The - power duration of 5.7 ms is equivalent to the spatial scale of ~250 
Debye lengths (A 0e ). The peak intensity E, of this wave packet is ~34.6 mVirT 


[Nicholson et al., 1978]. In the present case, this condition is 
easily satisfied, since the observed ^ of 8.5 x 10^ 4 , 2.6 x 
10 -3 , 3.2 x 10“ 3 and 1.3 x 10“ 3 are well above the esti- 
mated (k L X De ) 2 values of 1.9 x 1 0 4 , 2.2 x 10~ 4 , 1.2 x 
10~ 4 and 1.4 x 10 4 , for the events 1 to IV, respectively. 
This suggests that the strong turbulence processes can occur, 
and, therefore, in the context of present observations, we 
examine MI [Zakharov, 1972; Papadopoulos et al., 1974; 
Smith et al., 1979; Goldstein et al., 1979], and related soliton 
formation and Langmuir collapse [Zakharov, 1972; Nicholson 
et al., 1978; Goldman, 1984]. 

[ 26 ] For the events of the present study, the inequality 

<«> 


electrons by a factor of a = 1.3 to 1.7 [Alvarez et al., 1975 ;Lin 
et al., 1973], This implies that the beam speeds, corrected for 
the pitch angle scattering will be much more favorable for Ml 
and spatial collapse. 

[ 27 ] As far as the inequality ^ is concerned, it is also 
easily satisfied for all four events, since the observed are 
8.5 x 10“ 4 , 2.7 x 10“ 3 , 3.2 x 10“ 3 and 1.3 x"‘l0' 3 , 
whereas, — is ^5.5 x 1CV 4 . This indicates that the observed 

Mi 

wave packets are in the regime of supersonic modulational 
instability (MI) as shown schematically in Figure 18. The 
different regimes of Langmuir wave saturation mechanisms 
are discussed by Zakharov et al. [1985], 

[ 28 ] The growth rate of supersonic MI 


T_ = f WM 1/2 

LOpe \3 Ifli n e T e ) 


(5) 


is easily satisfied, since >,2 k D are of 2.3 x 1 0 3 m 1.7 x 

10~ 3 m _1 , 2.4 x 10~ 3 m~ and 3.8 x 1(V 3 m _I , whereas, the 
wave numbers of the beam excited Langmuir waves, k, are 
~1.3 x 10“ 3 m _1 , -1 x 10“ 3 m _1 , ~1.1 x 10“ 3 m _1 and 
~1.9 x 1 0" 3 m ~ 1 for the events 1 to IV, respectively (m e i are 
the electron and ion masses, respectively, and k D = A-). We 
have estimated the beam speeds for a smooth spiral. The pitch 
angle scattering is known to increase the path length of 


can be estimated as ^4 x 10 4 , ~7 x 10 4 , ~7.7 x 10 4 
and ~4.8 x 10~ 4 for the observed w ' f of 8.5 x 10 -4 , 2.7 x 
10 -3 , 3.2 x 10 -3 and 1.3 x 10 -3 corresponding to the 
events of Figures 4, 8, 12, and 16, respectively. On the other 
hand, the bandwidth of the beam-excited Langmuir waves 


A uj 

^pe 


3(£iAft,) 2 


M l 

h 


(6) 
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Figure 5. (a) The complete spectrum of the wave packet ofFigure 4 from 0 to 65 kHz, (b) the narrow spec- 
trum around f~f pe ~ 19.12 kHz, where L, D, and U correspond to the beam excited Langmuir wave, down- 
shifted sideband at ~ 1 9 kHz, and up-shifted sideband at ~ 1 9.3 kHz, respectively, and (c) the low frequency 
spectrum: the enhancement below 200 Hz corresponds to ion-sound waves. The sideband emissions are 
probably due to nonlinear interactions between the beam excited Langmuir waves and ion sound waves. 
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Figure 6. Dynamic spectrum of the second local type III radio burst (fast drifting emission from —5 MHz 
down to —20 kHz) and associated Langmuir waves (non-drifting emissions in the frequency interval 
10-14 kHz). 


can be estimated using the relation [Lin et al., 1986; Benz , 

2002] 

Aki_Anln2 

k L ~ v b IN' U 

Here A v b is the velocity dispersion of the beam, N — In p is 
the number of linear growth times before the onset of MI, and 
p is the ratio of the peak electric field amplitude to the thermal 
background. We estimate that p is — 1 0 4 for all four events. It 
is interesting to note that higher the N, the longer the time 
available for the linear beam-plasma instability to grow. For 
the events of the present study, we estimate the number of 
linear growth times N as >9 for p — 10 4 . Thus, for V — 9 and 
Av b - 0.1 v*, we obtain ^ — 3.9 x 10~ 3 . This value yields 
as -2.2 x 10“ 6 , -2.5 x 10 6 , —1.4 x 10' 6 and -1.6 x 
10 -6 for k L of — 1.3 x 10~ 3 m _1 , —1 x 10 -3 m _1 , —1.1 x 
10~ 3 m _1 and —1.9 x 10 -3 m -1 , and X De of 10.3 m, 14.0 m, 
9.8 m and 6.2 m, for the events I to IV, respectively. Thus the 
growth rates are much larger than the bandwidths of beam- 

^pe 

excited Langmuir wave packets This suggests that the 
pump waves are monochromatic enough for excitation of 
supersonic ML It is important to note that for the bandwidths 
of initial pump waves, one should not use the spectral widths 
from Figures 5b, 9b, 13b, and 17b since these are severely 
affected by the nonlinear effects. 

[ 29 ] The peak intensity also decides whether the wave 
packet is a collapsing Langmuir soliton of the type described 
by Zakharov [1972] and Nicholson et al. [1978] or not. For 


the wave packet to be a collapsing soliton, it should satisfy 
the condition [Thornhill and ter Haar, 1978; Gurnett et al., 
1981] 

> (A kX De f, (8) 

fie Te 

where A k = y is the wave number characteristic of the 
envelope. Assuming that the wave packet is convected in the 
solar wind, we can convert the measured timescales r of 
—5.7 ms, —7.9 ms, —10 ms and —7 ms into spatial scales S of 
— 25QX De , — 2&2X De — 5\2X De and — 355 using the rela- 
tion S — tv sw for v sw of 450 km s -1 , 500 km s -1 , 500 km s -1 
and 315 km s -1 , and X De of — 10.3 m, — 14.1 m, —9.8 m and 
—6.2 m, for the events I to IV, respectively. Thus, the 
observed of -8.5 x 10' 4 , -2.7 x 10“ 3 , -3.2 x 10“ 3 
and— 1.3 x 10 -3 are greater than (AkX De ) 2 of —6.3 x 10“ 4 , 
—5 x 10“ 4 , — 1.5 x 10“ 4 and— 3.1 x 10~ 4 obtained for the 
spatial scale S of — 250X De , — 282A £te , — 5\2X De and 
~?>55X De for events of Figures 4, 8, 12, and 16, respectively. 
This suggests that the observed wave packets are probably 
the Langmuir solitons collapsed to the spatial scales of a few 
hundred Debye lengths. 

3.1. Spectral Characteristics 

[ 30 ] As shown in Figure 1, MI, which is a four-wave 
interaction process excites a low frequency ion density per- 
turbation of frequency and wave number (Cl, q). This per- 
turbation can be either a freely propagating ion sound wave 
or a strongly damped quasi-mode. It can beat with two of the 
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Figure 7. Frequency-time spectrogram of the Langmuir waves associated with the type III burst of Feb- 
ruary 11, 2012. The Langmuir wave emissions, which are very bursty show substantial frequency spread- 
ing. The arrow shows the Langmuir wave burst corresponding to the current TDS event. 


initial pump waves with wave numbers, k L and frequencies 
f L ~ f pe (beam-excited Langmuir waves), and produce the 
down-shifted (Stokes) and up-shifted (anti-Stokes) side- 
bands with wave numbers and frequencies of ( k L — q ) and 
{f L — f 1), and ( k L + q ) and (f L + f 1), respectively. Thus, the 
Fourier spectrum resulting from such a four wave Ml should 
contain the main emission peak at f L ~ f pe together with two 
sidebands at (f L — fl) and {f L + f l), and a low frequency 
component. In the present case, the spectral peaks seen in 
Figures 5b, 9b, 13b, and 17b probably correspond to the 
down-shifted sideband (Stokes (D)), beam-excited Langmuir 
(L) and up-shifted sideband (anti-Stokes (U)) modes, 
respectively, and the low frequency spectral enhancements 
seen in Figures 5c, 9c, 13c, and 17c correspond to ion 
acoustic waves. This suggests that the observed spectral 
structures of the wave packets show one-to-one correspon- 
dence with the characteristic spectral signatures of the MI. 

[ 31 ] The requirement for the observed spectral compo- 
nents to be associated with MI is that they should satisfy the 
frequency, wave number and phase matching conditions: 


% = /d +fu 

(9) 

2 k L = k D + ku 

(10) 

2 4>l — ‘t’D + 0£/> 

(11) 


where, the subscripts L, D and U correspond to the beam- 
excited Langmuir, down-shifted sideband (Stokes) and up- 
shifted sideband (anti-Stokes) modes, respectively. 

[ 32 ] In order for the frequency matching condition to be 
satisfied, the down-shifted (Stokes) and up-shifted (anti- 
Stokes) sidebands should be symmetric with respect to the 
beam-excited Langmuir wave, and their frequency shifts 
Af = f — fpe should be in the frequency ranges of the 
respective low frequency enhancements. 

[ 33 ] For December 19, 2011 event (Figure 5b), the fre- 
quency shifts of the down-shifted and up-shifted sidebands 
are reasonably symmetric with respect to the Langmuir wave 
pump, being ~120 Hz and ~ 180 Hz, respectively. These 
frequency shifts are in good agreement with the observed 
frequencies of ion sound waves of <200 Hz (Figure 5c). For 
February 11, 2012 event (Figure 9b), the frequency shifts of 
the down-shifted and up-shifted sidebands with respect to 
Langmuir waves are ~220 Hz and ~250 Hz, respectively, 
which are in agreement with the observed frequencies of ion 
sound waves of ~250 Hz. For July 18, 2009 event 
(Figure 13), the frequency shifts of the D and U modes 
(Figure 13b) are ~250 Hz and ~220 Hz, which are again in 
good agreement with the observed frequencies <250 Hz of 
the low frequency enhancement as seen in Figure 13c. The U 
mode in this case is weaker than the D mode, probably 
because it suffers higher damping rate in comparison with 
that of the D mode. 
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Figure 8. One of the most intense Langmuir wave packets captured by the Time Domain Sampler (TDS) 
during the type III event of Figure 6. The A power duration of 7.9 ms of this event is equivalent to the spa- 
tial scale of 282A/} e . The peak intensity E L of this wave packet is ~44.7 mVm -1 . 


[34] In the case of TDS event of December 5, 2006, the 
spectral peak at ^31.6 kHz in Figure 17b corresponds to the 
beam-excited Langmuir wave, the spectral peaks below this 
frequency correspond to Stokes modes propagating anti- 
parallel to the beam, and the spectral peaks at frequencies 
higher than 31.6 kHz correspond to anti-Stokes modes 
propagating parallel to the beam. The connection between the 
structure in the spectrum of this event and the MI is as fol- 
lows. The MI excites low frequency ion density perturbations 
of frequency and wave number (Q, q), which can beat with 
two of the initial beam excited Langmuir waves, producing 
high frequency sidebands with wave numbers k L — q (Stokes 
mode) and Ay, + q (anti-Stokes mode). These Stokes and anti- 
Stokes modes can again couple to the pump wave and 
regenerate the initial low frequency density perturbation at <7; 
they can also couple to each other to generate another low 
frequency density perturbation at 2 q. This new density per- 
turbation can beat with the pump to produce another pair of 
Stokes and anti-Stokes waves at Ay ± 2 q which beat together 
to produce another density perturbation at 4 q, which, in turn, 
will also generate new Stokes and anti-Stokes waves. 
A hierarchy of low frequency density perturbations together 
with higher order sidebands are thus generated. This process 
of producing higher order sidebands can proceed as long as 
the level of ion density fluctuations are above the thermal 
level. The frequency shifts of all the Stokes and anti-Stokes 


modes are in the frequency range of the low frequency 
enhancement seen in Figure 17c. 

[35] We can estimate the wave numbers of the down- 
shifted and up-shifted sidebands and examine whether they 
satisfy the required matching condition as follows. The wave 
numbers k of the sidebands can be estimated using the 
expression for the Doppler shift: 

A/" = y-r*— (A'A Be )cos 0, (12) 

47 TADe 


where 8 is the angle between k and y7 w , i.e., 8 = 0 and 8 = tt 
correspond to the anti-Stokes and Stokes modes propagating 
away from and toward the Sun, respectively. The pump 
Langmuir waves, the sidebands and ion sound waves should 
satisfy the conservation relation 

k = k L ±q, (13) 


where the plus and minus signs correspond to anti-Stokes 
and Stokes modes, respectively. Since the wave numbers of 
the pump waves Ay are usually much smaller than q, the \q\ 

should be - 


[36] In the case of TDS event of Figure 4, the frequency 
shifts of the D and U modes of ~120 Hz and ~180 Hz as 
seen in Figure 5b are <200 Hz, which is the frequency range 
of the low frequency enhancement seen in Figure 5c. 
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Figure 9. (a) The complete spectrum of the wave packet of Figure 8 from 0 to 65 kHz, where the primary 
peak at ~14 kFLz corresponds to beam excited Langmuir wave, (b) the narrow logarithmic spectrum 
around / ~ f pe ~ 14 kHz, where the spectral peaks L, D, and U correspond to the beam excited Langmuir 
wave, down-shifted sideband ~ 13.78 kHz, and up-shifted sideband ~14.25 kHz, respectively, and (c) the 
low frequency linear spectrum: the enhancement below 250 Hz corresponds to ion-sound waves. The non- 
linear interactions between the beam-excited Langmuir waves and the ion sound waves observed are prob- 
ably responsible for the sideband emissions. 
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Figure 10. Dynamic spectrum of the third local type III radio burst (fast drifting emission from ~5 MHz 
down to ~20 kHz) and associated Langmuir waves (non-drifting emissions in the frequency interval 19— 
22 kHz). 
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Figure 11. Frequency-time spectrogram of the Langmuir waves associated with type III burst of July 1 8, 
2009. The arrow shows the Langmuir wave burst corresponding to the current TDS event. 
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Figure 12. One of the most intense Langmuir wave packets captured by the Time Domain Sampler 
(TDS) during the type 111 event of Figure 10. The 1 power duration of 10 ms of this event is equivalent 
to the spatial scale of 5\2\ De . The peak intensity E r of this event is ~70.9 mVnH 1 . 


Substituting these measured values A f ~ 120 Hz and A f= 
180 Hz in equation (12), we estimate kX De as 0.017 and 
~0.026, for the Stokes and anti-Stokes modes with 9 = 180° 
and 9 = 0° , respectively. These values indicate that the pump 
Langmuir waves with k L \ De ~ 0.014 are very efficiently 
converted into those of forward and backward propagating 
daughter Langmuir waves with large wave numbers. As far 
as the ion sound waves seen in Figure 5 c are concerned, 
since their phase velocities are usually well below v sw , their 
frequency spectrum in the spacecraft frame of reference is 
solely determined by the Doppler shifts. Thus, the upper 
limit of the wave numbers of these waves can be estimated 
as q = — 2.8 x 10~ 3 m~' and q\o e — 0.029 for S2 = 

200 Hz and v sw = 450 km s _1 . Thus, the wave numbers of 
the observed low frequency waves are comparable to those 
of the sideband emissions. 

[ 37 ] In the case of TDS event of Figure 8, the frequency 
shifts of the D and U modes ~220 Hz and ~250 Hz are in 
agreement with the observed frequencies of ion sound waves 
of ~250 Hz. Substituting these measured values A/ ~ 
220 Hz and A f= 250 Hz in equation (12), we can estimate 
k\ De as 0.039 and ^0.044, for the Stokes and anti-Stokes 
modes, respectively, and compare them with k L X De ~ 0.015. 
The upper limit of the wave numbers of the ion sound waves 
of Figure 9c can be estimated as q = — 3.1 x 10~ 3 nH 1 

and qX De — 0.044 for O = 250 Hz and v sw = 500 km s -1 . 


Thus, the wave numbers of the observed low frequency 
waves are comparable to those of the sideband emissions. 

[ 38 ] In the case of July 18, 2009 event (Figure 12), for the 

Stokes and anti-Stokes modes with 9 = 180° and 9 = 0°, we 
can plug the measured values of A f ~ 250 Hz and 220 Hz in 
equation (12) and estimate kX De as 0.031 and ~0.027, 
respectively. These values of kX De are much higher than 
k L X De ~ 0.015. As far as the ion sound waves seen in 
Figure 13c are concerned, the upper limit of their wave 
numbers are estimated as q = — 3.1 x 10~ 3 nH 1 and 

qXoe ~ 0.031 for n = 250 Hz and v sw = 500 km s -1 . Thus, 
the wave numbers of the observed low frequency waves are 
comparable to those of the sideband emissions, which indi- 
cates that the observed low frequency waves are probably 
involved in parametric interaction with the Langmuir waves. 

[ 39 ] In the case of Dec 5, 2006 event (Figure 16), we can 

estimate k L X De ~ 1.2 x 10~ 2 for Ay is ~1.9 x 10~ 3 m _1 and 
X De ~ 6.2 m. Plugging the maximum measured A/ of 
200 Hz in equation (12), we estimate that the wave numbers 
k of the anti-Stokes and Stokes modes should be as large as 
~2Ay. Thus, the spectral observations of Langmuir wave 
packets show that the beam excited Langmuir waves pre- 
sented in Figure 16 with k L X De of ^0.012 are very efficiently 
converted into those of the forward and backward propa- 
gating daughter waves with large wave numbers up to kX De 
~ 0.025. As far as the wave numbers of low frequency 
waves (Figure 17c) are concerned, we can estimate q = — 
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Figure 13. (a) The complete spectrum of the wave packet of Figure 12 from 0 to 65 kHz, where the dom- 
inant spectral peak at ~20. 1 kHz corresponds to the Langmuir wave excited by the electron beam at local 
electron plasma frequency, f pe , (b) the narrow spectrum around the main Langmuir wave peak at f ~f pe ~ 
20.1 kHz, where L, D, and U correspond to the beam excited Langmuir wave, down-shifted sideband at 
~ 19.85 kHz, and up-shifted sideband peak at ~20.32 kHz, respectively, and (c) the low frequency spec- 
trum: the enhancement below 250 Hz probably corresponds to ion-sound waves. 
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Figure 14. Dynamic spectrum of the fourth local type 111 radio burst (fast drifting emission from ~5 MHz 
down to ~32 kHz) and associated Langmuir waves (non-drifting emissions in the frequency interval 
28-32 kHz). 
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Figure 15. Frequency-time spectrogram of the Langmuir waves associated with the type 111 burst of 
December 5, 2006. The Langmuir wave emissions are very clumpy and show substantial frequency 
spreading. The arrow shows the location of the current TDS event. 
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Figure 16. One of the most intense waveforms of Langmuir waves observed by the Time Domain Sam- 
pler (TDS) during the type III event of December 5, 2006. The f power duration of 7 ms is equivalent to 
the spatial scale of 355\ De . The peak intensity E L of this wave packet is ^69.7 mVm -1 . 


4 x 10~ 3 m~' and q\o e — 0.025 for S2 = 200 Hz and v sw = 
315 km s -1 . These wave numbers of low frequency waves 
are comparable to those of sideband waves. 

[ 40 ] Thus, the sidebands and low frequency enhancements 
contained in the spectra of TDS events of Figures 4, 8, 12, 
and 16 probably correspond to the daughter waves excited 
by the Ml. 

3.2. Trispectral Analysis 

[ 41 ] The conventional spectral analysis does not yield 
complete information about the nonlinear processes from the 
waveform data. To extract the phase coherence amongst 
spectral components of the wave packets, the higher order 
spectral analysis techniques, namely, the bispectral and tris- 
pectral analysis techniques should be used. The bi-spectral 
analysis can extract the nonlinear characteristics between 
three spectral components with frequency sum and difference 
[Kim and Powers, 1979]. Since bi-spectral analysis is related 
to the skewness of a signal, it can detect only the asymmetric 
non-linearities, such as the electrostatic decay [ Henri et al. , 
2009; Balikhin et al, 2001; Walker et al., 2002, 2003] and 
harmonic generation [Bale et al., 1996], The amplitude of the 
bispectrum at the bifrequency ( F k , FI) measures the amount 
of coupling between the frequencies, F k , Fj and F k+t . The tri- 
spectral analysis, on the other hand, is used to demonstrate 
the occurrence of cubic nonlinear interactions in the system. 
The tricoherence spectra can isolate the individual cubic 


interactions between quartets of spectral components, espe- 
cially those associated with oscillations at the power spectral 
primary peak frequency and its sidebands, as well as their 
combination tones. 

[ 42 ] A cumulant based trispectral method has been 
developed and applied to synthetic [Kravtchenko-Berejnoi et 
al., 1995a, 1995b; Lefeuvre et al., 1995] and simulated data 
[Soucek et al., 2003]. The expression for the trispectrum can 
be written as [Kravtchenko-Berejnoi et al., 1995a] 

T(k,l,p) = E[x k X,X' p xi\ - N(k,l,p,q ), (14) 

where X k , X h X p and X q are the complex Fourier components 
of the signal at frequencies F k , F h F p and F q = F k + F/ — F p , 
and N(k, l, p, q) = E[X k X,\E[XX] + E[X k 4]E[XX,] + E 
[Ay^]7?[A/A^]. Here E{\ is the expectation operator. In the 
present case, the subscripts k and / refer to beam excited 
Langmuir waves, and p and q refer to down-shifted (Stokes) 
and up-shifted (anti-Stokes)sidebands, respectively. 

[ 43 ] The normalized trispectrum is usually referred to as 
the tricoherence. It is more useful because it eliminates the 
dependence of trispectrum on the amplitudes of the signals. 
The tricoherence quantifies the phase coherence amongst the 
spectral components of the wave packet, i.e., it quantifies the 
fraction of the total product of powers at the frequency 
quartet, (F k , F/, F p , F q = F k + F p — Fj), that is owing to 
cubicly phase-coupled modes. The expression for the square 
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Figure 17. (a) The complete spectrum of the wave packet of Figure 16 from 0 to 65 kHz, where the dom- 
inant spectral peak at ^31.6 kHz corresponds to the beam excited Langmuir wave, (b) the narrow spec- 
trum around / ~ f pe ~ 31.6 kHz, where the main peak is denoted as L, the down-shifted (Stokes) 
sidebands are denoted as D\ and D 2 , and up-shifted (anti-Stokes) sidebands are denoted as U\, U 2 and 
U 2 . (c) The low frequency spectrum, which shows a clear enhancement below 200 Hz corresponding to 
ion-sound waves. 
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Table 1. Relevant Parameters 


Parameters 

19/12/2011 

11/2/2012 

18/7/2009 

5/12/2006 

Electron Density, n e m -3 

4.5 x 10 6 

2.4 x 10 6 

5 x 10 6 

1.2 x 10 7 

Solar wind velocity, v sw km s -1 

450 

500 

500 

315 

Electron Temperature, T e ( K) 

10 5 

10 5 

10 5 

10 5 

Debye length, \ De (m) 

10.3 

14 

9.8 

6.2 

Electron plasma frequency,^ (kHz) 

19.12 

14 

20.1 

31.6 

Electron Thermal velocity, vy e (ms -1 ) 

1.23 x 10 s 

1.23 x 10 6 

1.23 x 10 6 

1.23 x 10 6 

Beam velocity, v h 

0.3c 

0.28c 

0.37c 

0.35c 

Beam width, ^ 

~0. 1-0.2 

-0. 1-0.2 

~0. 1-0.2 

~0. 1-0.2 

Beam resonant wave number, (m -1 ) 

1.3 x 10~ 3 

1 x 10~ 3 

1.1 x 10~ 3 

1.9 x 10~ 3 

Peak wave amplitude, E L (mVm -1 ) 

34.6 

44.7 

70.9 

69.7 

Normalized energy density, 

8.5 x 10~ 4 

2.6 x 10~ 3 

3.2 x 10~ 3 

1.3 x 10~ 3 

—level duration (ms) 

5.7 

7.9 

10 

7 

l-level size (Ace) 

250 

282 

512 

355 

Peak tricoherence 

0.35 

0.58 

0.32 

0.32 


tricoherence can be written as \Kravtchenko-Berejnoi et al . , 
1995a]: 




\T{k,l,p) | 2 

(E[\x k x,x;x;\]+\N(k,i,p,q)\y 


(15) 


method involves the division of the data record into M seg- 
ments; an appropriate window is applied to each segment to 
reduce leakage; the trispectrum as well as tricoherence are 
computed for each segment by using the DFT; finally, the 
trispectrum and tricoherence are averaged across segments to 
reduce the variance of the estimator. 


A unit value of the tricoherence indicates perfect coupling, 
a zero value indicates no coupling, and any value between 
zero and one indicates partial coupling. 

[ 44 ] The tricoherence is zero for a Gaussian process. Even 
for a Gaussian process, due to statistical fluctuations, the 
estimate of the tricoherence from a finite data record will not 
be zero. Therefore, the method of periodograms is usually 
used to estimate the trispectrum and tricoherence. This 



I - Decay 

II - Modulational Instability 

III - Subsonic Modulational 

IV - Supersonic Modulational 

V - Modified Decay 

Figure 18. Different regimes of Langmuir turbulence. 


STEREO A - 19-Dec-2011: 13:53:46.131 





0.0 0.1 0.2 0.3 0.4 

t 2 

Figure 19. The three dimensional representation of the tri- 
coherence spectrum t z (F k , F h F p ) of the TDS event of 201 1 
December 19 (Figure 4). The tricoherence r exhibits peak 
value of ^0.35 at (~19.25, ~19, ~18.75) kFIz, which quan- 
tifies the phase relation 2 <f> L = 4 >d + where <j> L , (f) D and <j>u 

are the phases of the beam-excited Langmuir wave at 
~ 19. 12 kHz, Stokes (~19 kHz) and anti-Stokes 
(~19.3 kHz) modes, respectively. 
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STEREO B - ll-Feb-2012:20:24:41.650 



0.0 0.1 0.2 0.3 0.4 0.5 

t 2 

Figure 20. The 3-D representation of the tricoherence 
spectrum t 2 (F k , F b F p ) of the TDS event of Figure 8 (2012 
February 11). The peak tricoherence of —0.58 at frequency 
quartet ( F k — 14.5, F/ — 14.25, F p — 13.75, F q = F k + F p — 
F[ — 14.125) kFlz quantifies the phase coherence between 
two beam excited Langmuir waves and, a down-shifted and 
an up-shifted sidebands. 

[ 45 ] As seen from equation (14), the trispectrum estimator 
is symmetric with respect to permutations of its arguments 
F k ,F l and F p . The principal domain for the interaction of the 
type F k + Fi ~ F p + F q is determined as [Kravtchenko- 
Berejnoi et al., 1995a] 0 < F k < f N , 0 < F t < F k , 0 < F p < 
F h and F p < F k + F/ — F p <f N , where f N is the Nyquist 
frequency. 

[46] Thejappa et al. [2012a] reported the observations of a 
highly localized Langmuir wave packet captured by the TDS 
during one of the local type III events. This wave packet was 
characterized by a short duration of —3 ms, peak intensity 
exceeding the threshold for the OTS1 as well as spatial col- 
lapse, and the spectrum, containing the signatures of OTSI, 
namely, an intense peak corresponding to beam-resonant 
Langmuir wave, two sidebands corresponding to down- 
shifted and up-shifled daughter Langmuir waves, and a low 
frequency enhancement corresponding to daughter ion sound 
waves; the frequencies and wave numbers of these spectral 
components satisfied the resonance conditions of OTSI. 
Since the spectral methods are phase blind, Thejappa et al. 
[2012b] had developed cumulant based trispectral analysis 
tools to extract the phase information. Using those tools, 
these authors have clearly demonstrated that the spectral 
components of the observed wave packet are coupled to each 
other with a high degree of phase coherency as expected of 
OTSI. 


[ 47 ] In this study, we will use the trispectral analysis tools 
developed by Thejappa et al. [2012b] to extract the phase 
information amongst the spectral components of the wave 
packets of Figures 4, 8, 12, and 16. For this purpose, as 
expected of the periodogram method, we divide the time 
series into M = 16 segments with segment length N = 1000 
(0.004 s). We apply the Flamming window for each segment, 
and calculate the tricoherence spectrum as a function of three 
frequencies F k , F h and F p . 

[48] The true values of tricoherence are known to be zero 
for Gaussian noise \Hasselman et al., 1963; Nikias and 
Raghuveer, 1987; Dalle Mode and Hinich, 1989]. The 
90% significance level for estimate of tricoherence for 
Gaussian noise is given by [Chandran et al., 1994] 

= ( 16 ) 

In the present case, for M = 16, we obtain t 2 c — 0.144. 
Therefore, the tricoherence r values above 0.144 are sig- 
nificant at the 90% significance level for M blocks of the 
data, i.e., for 2 M degrees of freedom. Although any trico- 
herence value between zero and one indicates a partial 
coupling of four waves, we consider only those values which 
are greater than 0.144 as significant. 


STEREO A - 18-Jul-2009:02:42: 17.159 



0.00 0.05 0.10 0.15 0.20 0.25 0.30 


t 2 

Figure 21. The 3-D tricoherence spectrum t 2 (F k , F b F p ) of 
the TDS event of Figure 12 (2009 July 18). The peak tricoher- 
ence of -0.32 at (F k - 20.25 F, - 20, F p - 19.75, F q = 
F k + F p F/ — 20) kHz, quantifies the phase relation 
2 (j) L = (j)D + <t>u, where 4>l, <1>d and cj) V are the phases of the 
beam-excited Langmuir, Stokes and anti-Stokes modes, 
respectively. 
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STEREO B - 05-Dec-2006: 12:09:55.526 




t 2 


Figure 22. The 3-D tricoherence spectrum f(F k , F,, F p ) of 
the TDS event of Figure 16 (2006 December 5). The peak 
coherence —0.32 at the frequency quartet (F k — 31.5, F/ — 
3 1 .5, F p — 3 1 .25, F q = F k + F t — F p — 3 1 .75) kHz quantifies 
the phase relation 2 <j> L = (f> D + fy, where (j) L , <j> D and fy are 
the phases of the beam-excited Langmuir, one of the Stokes 
and one of the anti-Stokes modes, respectively. 


[ 49 ] Since it is difficult to visualize the tricoherence results 
in a 3-D space, usually one frequency is fixed and the trico- 
herence spectrum is displayed as a function of two frequen- 
cies. However, in the present study, we display the results in 
3-D space, i.e., we display the tricoherence t 2 spectrum as a 
function of three frequencies, F k , F t and F p . We restrict the 
frequencies to narrow relevant frequency ranges, centered 
around the frequency of the primary spectral peak, which in 
the present case is the local electron plasma frequency, f pe . 
We construct the 3-D spectrum t 2 (~f D < F k < ~fy, ~f D <F/< 
~fu, ~fo < F p < ~fy), where f D and fu correspond to the 
frequencies of down-shifted and up-shifted sidebands, 
respectively. Since the frequency resolution in the present 
case is —250 Hz, the considered frequencies may not exactly 
coincide with the frequencies of the observed spectral peaks. 

[ 50 ] In Figure 19, we present the 3-D tricoherence spec- 
trum of the TDS event of 2011 December 19, constructed 
using the tricoherence r values computed in the frequency 
ranges [19 kHz <F/< 19.75 kHz, 19 kHz < F k < 19.5 kHz, 
18.75 kHz <F p < 19.25 kHz]. As seen from this 3-D repre- 
sentation, the tricoherence exhibits a spectral feature of finite 
volume in 3-D frequency space. However, The tricoherence 
attains apeak value of —0.35 at (—19.25, — 19, — 18.75) kHz. 
This value of tricoherence, which is a statistically significant 
value by being greater than iy — 0.144 quantifies the phase 
relation 2 <j> L = <f) D + <f>y, where <p L , and fy are the phases 


of the beam-excited Langmuir wave at — 19. 12 kHz, Stokes 
(—19 kHz) and anti-Stokes (~ 19.3 kHz) modes, respectively. 

[ 51 ] For 2012 February 11 event, we present the 3-D tri- 
coherence spectrum in Figure 20. This 3-D spectrum is 
constructed using the computed f values in the frequency 
range [14.25 kHz<.F*< 14.75 kHz, 14kHz<F,< 14.25 kHz, 
13.75 kHz < F p < 14.25 kHz]. Although, the tricoherence 
spectral feature has an extended shape in 3-D frequency 
space, it exhibits a peak of —0.58 at (—14.5, — 14.25, 
— 13.75) kHz. This peak quantifies the phase coherence 
between two beam excited Langmuir waves at —14 kHz and 
a down-shifted sideband at — 13.75 kHz and an up-shifted 
sideband at — 14.25 kHz. 

[ 52 ] For the 2009 July 18 event, the 3-D tricoherence 
spectrum is presented in Figure 21. The tricoherences used 
in this figure are computed in the frequency ranges [20 kHz 

< F k < 20.75 kHz, 20 kHz <F/< 20.25 kHz, 19.75 kHz < F p 

< 20.25 kHz]. In this case, the t 2 attains a peak of —0.32 at 
(—20.25, —20, — 19.75) kHz. Again, since it is statistically 
significant, it quantifies the phase relation 26/ = 4>d + (f>y ■ 

[ 53 ] In Figure 22, we present the 3-D tricoherence t 2 
spectrum of the 2006 December 5 event. In this case, the 
tricoherence r attains a peak of —0.32. This is computed in 
the frequency range [31.625 kHz < F k < 31.875 kHz, 
31.5 kHz <F/< 31.875 kHz, 31.25 kHz < F p < 3 1 .75 kHz], 
The corresponding frequency quartet is (^3 1.5, —31.5, 
—31.25, —31.75) kHz. This tricoherence quantifies the 
phase relation 2 <f> L — (p D + <f>y. 

[ 54 ] The results of tricoherence estimates can also be dis- 
played as 2-D cross-sections at different frequencies. For 
example, for 2006 December 5 event, we present the cross 
sections of the tricoherence spectrum in Figure 23. The top 
panel shows the cross-section at F p =f D = 3 1 .25 kHz, which 
shows a clear tricoherence spectral peak at (^3 1.5, —31.5, 
—31.25) kHz. The peak value of t in this case is —0.32. 
This value agrees with that of 3-D analysis. This peak tri- 
coherence quantifies the phase relation 2(f) L = (f>o + <f>u, 
where (f> k , (f>o and (f> v are the phases of the beam-excited 
Langmuir wave (~31.5 kHz), one of the Stokes 
(~31.25 kHz) and one of the anti-Stokes (~3 1 .75 kHz) 
modes, respectively, i.e., the frequency quartet in this case is 
(~31.5, ~31.5, ^3 1 .25, ~31.75) kHz. The maximum tri- 
coherence I 2 ~ 0.32 in this case is statistically significant. 

[ 55 ] In the second panel, we present the cross-section of the 
tricoherence spectrum at F p =f L = 31.50 kHz. It shows two 
comparatively weaker spectral features, with peak values of 
-0.21 at (-31.75 kHz, -31.625 kHz, -31.50 kHz) as well 
as at (—32, —31.875, —31.50) kHz. These relatively smaller 
tricoherence values reflect the weaker four wave interactions 
involving the sidebands themselves. In the bottom panel, we 
present the tricoherence cross-section at F p =fy= 30.75 kHz, 
which also shows a weak peak of —0.21 at (—31.75, 
—31.875, —32.125) kHz. This again reflects the four wave 
interactions involving the sidebands themselves. However, 
the tricoherence spectral peak t — 0.32 (top panel) is more 
significant in the statistical sense. Therefore, the dominant 
four wave interaction in this case is 2 '.f L —*f D + fu- One should 
note that, the structure in the spectrum of this TDS event as 
seen in Figure 17b with multiple sidebands and the trico- 
herence peaks in the second and third panels indicate that the 
four wave interactions between the sidebands also can be 
possible. 
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Figure 23. The cross-sections of the tricoherence spectrum of the TDS event of Figure 16, at (top) 
F p =f D = 3 1 .25 kHz, (middle) F p =f L = 31.50 kHz, and (bottom) F p =fu= 3 1 .75 kHz, where fo,fL an d fu ar e 
the frequencies of the down-shifted sideband, beam-excited Langmuir and up-shifted sideband modes, 
respectively. 
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[56] Thus, the significant tricoherences corresponding to 
2 <f> L = 4 >d + fit; in the tricoherence spectra at Stokes, anti- 
Stokes and beam-excited Langmuir wave frequencies provide 
evidence for the modulational instability type of four-wave 
interactions in all the four events presented in this study. 

4. Conclusions 

[57] We have presented the observational evidence for the 
strong turbulence processes, namely, the supersonic mod- 
ulational instability (MI), which is also known as the oscil- 
lating two stream instability (OTSI) and collapsing 
Langmuir solitons in the source regions of solar type III 
radio bursts. The high time resolution observations from the 
STEREO WAVES experiment show that the Langmuir 
waves excited by the electron beams occur as localized wave 
packets with durations less than ~10 ms and peak intensities 
exceeding the strong turbulence thresholds. These localized 
Langmuir wave packets satisfy the criterion of solitons col- 
lapsed to spatial scales of a few hundred Debye lengths. The 
spectra of these wave packets consist of intense peaks 
corresponding to beam-resonant Langmuir waves, two or 
more sidebands corresponding to Stokes and anti-Stokes 
modes, and low frequency enhancements below a few hun- 
dred Hz corresponding to ion sound waves. The frequencies 
and wave numbers of these spectral components satisfy the 
resonance conditions of the modulational instability. The 
tricoherence values are computed using the trispectral anal- 
ysis techniques, which show that the spectral components of 
the wave packets are coupled to each other with a high degree 
of coherency, as expected of the ML The observed high 
intensities, short scale lengths, sideband spectral structures 
and low frequency enhancements of the wave packets, and 
high levels of tricoherence amongst their spectral compo- 
nents strongly suggest that the wave packets are collapsing 
Langmuir solitons, and they are in the supersonic modula- 
tional instability regime. The implications of the observations 
presented in this study are that ( 1 ) the strong Langmuir tur- 
bulence processes occur in the source regions of solar type III 
radio bursts, (2) the supersonic modulational instability (MI) 
and spatial collapse probably play important roles in beam 
stabilization as well as conversion of Langmuir waves into 
escaping radiation at the fundamental and second harmonic 
of the electron plasma frequency, and (3) the Langmuir col- 
lapse in the case of type III bursts probably follows the 
supersonic modulational instability as suggested by 
Zakharov [1972]. 
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